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SUMMARY 

A method is developed for the minimization of gear tooth surface deviations 
between theoretical and real surfaces to improve the precision of surface manufac- 
ture. Coordinate measurement machinery is used to determine a grid of surface 
coordinates. Theoretical calculations are made for the grid points. A least-square 
method is used to minimize the deviations between real and theoretical surfaces by 
altering the manufacturing machine-tool settings. An example is given for a hypoid 
gear . 


INTRODUCTION 

The Gleason Works were pioneers in the application of coordinate measurements to 
improve the precision manufacturing of hypoid and spiral bevel gears (ref. 1). In 
aerospace applications, duplication of flight-qualified master gears is very impor- 
tant, and coordinate measurement has now become part of the normal production pro- 
cess. Methods to enhance and extend the use of this machinery can be very valuable 
to aerospace gear manufacturers. 

The approach developed in this paper enables one do determine deviations of a 
real surface from the known theoretical surface. This is accomplished by using 
coordinate measurements and minimizing the deviations to correct the previously 
applied machine-tool settings. The surface deviations are represented in the direc- 
tion of the normal to the theoretical surface. The coordinate measurements are 
performed by a machine with 4 or 5 degrees of freedom. In the case of 4 degrees of 
freedom, the probe performs three translational motions (fig. 1); the fourth motion, 
rotation, is performed by turning the table with the workpiece. The axis of 
rotational motion coincides with the axis of the workpiece. In the case of a 
5-degree-of-f reedom machine, the fifth degree of freedom is used to provide the probe 
deflections in the direction of the normal to the theoretical surface. The probe is 
provided with a changeable spherical surface whose diameter can be chosen from a wide 
range . 

The motions of the probe and the workpiece for coordinate measurements are com- 
puter controlled and for this purpose a grid, the set of points on the surface to be 
measured, must be chosen. There is a reference point, one point on the grid, that is 
necessary for the initial installments of the probe. There are two orientations of 
the probe installment that are applied for measurements of a gear (fig. 1(a)) and a 
pinion (fig. 1(b)), depending on the angle of the pitch cone. 
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Figure 1.— Surface measurement orientation for (a) gear 
and (b) pinion. 

The mathematical aspects of coordinate measurements will now be described 
(ref. 2): First, it is necessary to derive the equations for the theoretical sur- 

face. In many cases, this surface can be derived as the envelope to the family of 
g enera ting surfaces, namely the tool surfaces. Next, the results of coordinate 
measurements must be transformed into deviations of the real surface represented in 
the direction of the surface normal. Here, the surface variations are represented in 
terms of the corrections to the machine-tool settings. The surface deviations 
obtained from coordinate measurements and the surface variations determined by the 
corrections of machine-tool settings can be represented by an overdetermined system 
of linear equations. The number k of these equations is equal to the number of 
grid points, and the number of unknowns m is equal to the number of corrections of 
machine-tool settings (m « k) . The optimal solution to such a system of linear 
equations enables one to determine the sought-for corrections of machine-tool 
settings . 


EQUATIONS OF THEORETICAL TOOTH SURFACE E 

Considering that the theoretical surface can be determined directly, we repre- 
sent it in coordinate system S in two-parametric form as 

r t (u,0), n t (u,0) i 1 ) 

where r and n are the position vector and unit normal to the surface, respec- 
tively, and (u ,6) are the Gaussian coordinates (surface coordinates). 
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For the case when surface E t is the envelope to the family of generating 
surface E , we represent in S surface E and the unit normal n to E as 

c c t t t t 

(ref. 3) 

r t - [M tc ] r c(u c ,0 c ),f(u c ,0 c ^) =0 ( 2 ) 

n t * [L tc ]n c (u c ,0 c >'f<« c ' 0 c'^ “0 < 3 > 

where (u ,6 ) are the Gaussian coordinates of the generating surface E , and ^ is 
the generalized parameter of motion in the process for generation. The equation of 
meshing is 


f(u C '0.<p) = N 


(c) 


,(ct) 


(4) 


where N (c) is the normal to E , and v (ct) is the relative motion for a point of 
contact of E and E . Matrices [M ] and [L 1 describe the coordinate 

c t 1 tc J 4x4 1 tc J 3x3 

transformation from to S for a position vector and surface normal, respec- 

tively. Position vectors in three-dimensional space are represented with homogeneous 
coordinates . 


COORDINATE SYSTEMS USED FOR COORDINATE MEASUREMENTS 

Coordinate systems S^ and S t are rigidly connected to the coordinate measur- 
ing machine (CMM) and the workpiece being measured, respectively (fig. 2). The 



Figure 2. — Relationship between theoretical and measurement 
coordinate systems, (p, radial distance to point from axis of 
rotation.) 
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Figure 3.— Measurement grid on tooth surface. 


backface of the gear is installed flush with the base plane of CMM. The distance 1 
between the origins and is assumed to be known, but the parameter 6 of 

orientation must be determined (see the following section). The coordinate transfor- 
mation from S t to is represented by the matrix equation 

r n, = < 5 > 


MEASUREMENT GRID AND ESTABLISHMENT OF THE REFERENCE POINT 

The grid is a set of points on chosen as points of contact between the 

probe and E t (fig. 3). Fixing the value of z for the point of the grid and the 
value of, say, y (or x t ) , we can obtain the following equations: 

* hi* = ^i (i = 3/ *«>rk) (6) 

where k is the number of grid points. 

We consider h and 1 as given and solve equations (6) for ( u i /® i )* Then we 
can determine the position vectors and the unit normals for k points of the grid 
using the equations 

= z t ( u A> ] T ' (i = 1, . . . ,k) (7) 


nj. 11 = ln xt { u i' 0 i) n yt (u i ,0 i )n xt( u i' 0 iU T » ( i = 1 , . . . , k) 


The position vector for the center of the probe, if the deviations are zero, is 
represented by 


4 


( i) (i) 


(i) 


+ P n t 


( i = 1, . . . ,k) 


(9) 


where p is the radius of the probe sphere. 

The reference point 

(o) - - (o) Mo) . . (o) MO) , 






( 10 ) 


is usually chosen as the mean point of the grid. 

The center of the probe that corresponds to the reference point on is 

determined from equation (9) as 


(o) 


[x t (u (o, ,0 (o) )y t (u' O) r 0 O) )z t (u'°',0 lO) )] 


(O) flO), 


n(o) x , T 


( 11 ) 


where (u (o \0 (o} ) are known values. 


The coordinates of the reference center of the probe are represented in coordi- 
nate system S of the measuring machine by the matrix equation 


(o) c (o) 

\ =[M mt (5)]R t 


( 12 ) 


Equation (12) yields 


(O) (o 


*r<5, 


r';’ 


(°) (o) _ (o) -t Q ) 

z m (5, u< \e { 


(13) 


_ (O) (o) (o) 

Three equations (13) contain four unknowns: o, x m , y m , . To solve these equa- 

tions, we may consider that one of the coordinates of the reference point of the 

(O) 

probe center, say, y m , may be chosen equal to zero. This is accomplished by 
requiring the reference point to lie in the x^ - plane. The orientation of 

angle 6 is now established to satisfy this requirement, and all measurements are 
referenced from this location. Then equation system (13) allows one to determine 


r , (°) , . J , (o) (o) (o) 

0, x m and z (ref. 2). Coordinates x ,y * 0, z 
m m ' m J ra f m 

initial installment of the center of the probe. 


are necessary for the 
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MEASUREMENT OF THE DEVIATIONS OF THE REAL SURFACE 


The deviations of the real surface are caused by manufacturing errors, heat 
treatment, etc. Vector positions of the center of the probe for the theoretical 
surface and the real surface can be represented as follows: 

R m - r ( u , # ) + On (u f 9) 

R* = r(u,0) + Xn (u ,6) 


where r and n are the position vector and the unit normal to the theoretical 
surface and are represented in coordinate system S^ of the measuring machine; X 
determines the real location of the probe center and is considered along the normal 

* 

to the theoretical surface; R and R represent in S the position vector of the 

mm m 

probe center for the theoretical and real surfaces, respectively. Equations (14) 
and (15) yield 

R* - R = (X - p)n = Ann <16) 

m m ' I'm m 


and 


An = (R - R ) 


(17) 


The position vector R* is determined by coordinate measurements for points of the 
grid. Equation (17) determines numerically the function 

hn i = An i (u i , ^ i ) (i = 1, . . . ,k) (18) 

that represents the deviations of the real surface for each point of the grid. 


MACHINE TOOL SETTINGS TO MINIMIZE DEVIATIONS 

The procedure used to minimize the deviations can be represented in two stages: 
(1) determination of variations of theoretical surface caused by changes of applied 
machine-tool settings, and (2) minimization of deviations of real surface by appro- 
priate correction of machine-tool settings. 

We consider that the theoretical surface is represented in S t as 

r t = r t ( u »®*dj) ( j = 1, . . . ,m) (19) 

where parameters d^ are the machine-tool settings. The surface variation is 
represented by 


5r 


t 


9 r t x 9r t 

~&I 6u + 


m 5r\ 

« + E 

j-i 9d j 


( 20 ) 
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We multiply both sides of equation (20) by the surface unit normal n t and take into 
account that dr /9 9 • n t = 9r t /9u • n t = 0, since 9r t /00 and 0r t /9u lie in the 
plane that is tangent to the surface. Then we obtain 


5r t ‘ n t = E 

j = l 


9d, 


n* 


Sdj = £ a5d j 
j=i 


( 21 ) 


We can now consider a system of k linear equations in m unknowns (m « k) of the 
following structure: 


a il 5d l + a i2 5d 2 + • • * + a i m 5d n, = b l 


a kl 5d l + a K2 5d 2 + 


+ a. 5d m = b. 

Km m K. 


( 22 ) 


Here 


b i “ An i " (Ki - R n,i> * n n,i (23> 

where i designates the number of grid points; a . (s = l,...,k; j = 

represents the dot product of partial derivatives 9r t /9d^ and unit normal n t - The 
system (22) of linear equations is overdetermined since m << k. The essence of the 
procedure for minimization of deviations is determining unknowns 6d^ (j = l,...,m) 
that will minimize the difference between the left and right sides of equations (22). 
The solution employed the least-square method. The subroutine DLSQRR of IMSL MATH/ 
LIBRARY (ref. 4) was used to computerize the procedure. 


APPLICATION OF METHOD TO THE INSPECTION OF FORMATE HYPOID GEAR 

Each tooth side of a formate face-hobbed gear is generated by a cone, and the 
gear tooth surface is the surface of the generating cone. Two cones that are shown 
in figure 4(a) represent both sides of the gear space. The following equations 
represent in coordinate system S gear surfaces for both sides and the unit normal 
to such surfaces (fig. 4(b)): 


r 


c 


(r - s 0 sin a G )sin d Q 
(r - s G sin a G )cos 0 Q 


1 


(24) 
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sin a G 

- cos a G sin 9 q 

- COS a G COS 9 g 

where, r c is the position vector and n c the unit normal; r is the cutter tip 
radius; a is the cutter blade angle (a >0 for the concave side and a < 0 for 

G G G 

the convex side) . 

Figure 5 shows the installment of the generating cone on the cutting machine. 
Coordinate systems and are rigidly connected to the cutting machine and the 

gear being generated, respectively. Systems S^, S^, and S t are rigidly connected 
to each other since the gear is formate cut. To represent in S t the theoretical 
gear tooth surface E t and the unit normal to E t , we use the following matrix 
equations : 

r t ( S G'^G' d j ) = r M tc] r c< S G' 0 G) < 26 > 

n t< S G' 0 G' d j> = [ L tJ n c( S G' 0 G> 

where 

[M ] = [M ] [M ] 

1 tc J 1 to J 1 oc J 

cos 7 0 -sin 7 0 1 0 0 0 

0 1 0 0 0 1 0 -v 2 

sin 7 0 cos 7 -Ax 0 0 1 H, 

* m i m m 2 

0 0 0 1 J [o 0 0 1 

The surface Gaussian coordinates are s and 9 , and d,(7 , V , H , and Ax ) are 

G G j ' 'm 2 2 nr 

the machine-tool settings. 

The numerical example presented in this paper is based on the experiment that 
has been performed at the Dana Corporation (Fort Wayne, IN, U.S.A.). The initial 
deviations An for each side of the real tooth surface have been obtained by 
measurements on a coordinate measuring machine (fig. 6). The grid' for the meas- 
urement is formed by nine sections along the tooth length, each section having five 
points. The number of grid points k is therefore 45, and the reference point is at 
the middle of the grid, i.e., the third point of the fifth section. In the measure- 

(o) 

ment, the coordinate y of the reference point is chosen to be zero and the 
c m 

alignment angle 0 is determined from solving equation system (13). 

The minimization of deviations was performed in accordance with the algorithm 
described in MACHINE TOOL SETINGS TO MINIMIZE DEVIATIONS, and the results are 
illustrated in figure 7 and table I. 


(27) 


(28) 



( 25 ) 
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Figure 7. — Minimized deviations after corrections made to machine-tool settings. 


TABLE I. - RESULTS OF MINIMIZATION 

[Pressure angle, a Q = 21.25°; cutter diameters « 9 in.; 
point width of cutters = 0.08 in.] 


Machine-tool 

settings 

Machine-tool setting parameters 

V 

mm 

H 2' 

mm 

rad 

Ax , 

rn 

mm 

Initial 

Corrected 


27.466600 

27.21603 

1.059816 

1.06437 

0.009677 

-0.53343 


CONCLUSION 

A general approach for a computerized determination of deviations of a real sur- 
face from the theoretical one based on coordinate measurements has been proposed. An 
algorithm for computerized minimization of deviations by corrections of initially 
applied machine-tool settings has been developed. The approach is illustrated with 
the example of the tooth surface of a hypoid formate gear. 
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